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摘要 : 导航 信号 作为 射电 干扰 信号 影响 工 波 段 射电 观测 时 ， 与 雷达 、 通 信和 领域 的 干扰 
不 同 ， 其 功率 通常 小 于 噪声 功率 ， 使 得 利用 干扰 信号 相关 性 的 传统 抗 干扰 方法 效果 不 佳 ， 且 
缺乏 有 效 的 评价 指标 。 针 对 以 上 问题 ， 依 据 射 电信 号 处 理 流程 构建 了 射电 抗 干扰 系统 框架 ， 
= 并 从 射电 观测 数据 应 用 出 发 ， 给 出 了 干扰 抑制 度 和 信号 损失 度 两 个 指标 。 同 时 提出 在 精确 补 
偿 时 延 后 利用 抑制 干扰 子 空间 骂 声 分 量 的 自 适 应 这 办 相 消 方法 来 消除 导航 信号 的 影响 ;仿真 
N 结果 表明 ， 本 方法 能 有 效 抑 制导 航 信 号 对 射电 观测 的 影响 ， 且 性 能 优 于 传统 基于 特征 子 空 间 
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随 着 射电 天 文 的 发 展 ， 国 际 通讯 联盟 所 划分 的 射电 观测 频段 已 不 能 满足 射电 天 文 研究 的 需求 ， 观 
测 保护 频段 以 外 的 射电 源 ， 常 常 受到 干扰 信号 的 影响 ， 极 大 限制 了 射电 源 的 观测 质量 。 导 航 信 号 是 工 
频段 射电 观测 的 常见 干扰 ， 根 据 全 球 定位 系统 和 北斗 导航 卫星 的 ICD 文件 可 知 " 2 ， 全 球 定位 系统 导航 
> 音 号 的 主要 频率 是 L1: 1575. 42 + 1.023 MHz, 12; 1 227.6 + 10.23 MHz 等 ， 北 斗 导 航 信号 的 主要 频率 
-一 是 : Bl: 1556.098 42.046 MHz, B3: 1268.52+10.23 MHz 等 ,而 这 些 频 率 与 红 移 的 分 子 谱 线 、 脉 冲 
e 星 等 射电 源 的 观测 频段 重合 。 同 时 ， 根 据 ICD 文件 中 导航 信号 的 最 小 电 平 可 以 换算 其 功率 通 量 谱 密 度 为 
c -196~ -201 dBw/(m?/Hz) ， 而 在 同样 带宽 下 1 mJy 的 射电 信号 功率 通 量 谱 密度 为 -290 dBw/(m?/Hz) , 
因此 在 接收 时 导航 信号 比 射 电信 号 强 89~91 dB。 此 外 ， 天 线 接收 系统 对 接收 信号 也 有 影响 ， 由 于 导航 信 
号 往往 从 旁 狼 进 入 射电 接收 机 ， 经 过 放大 、 滤 波 等 过 程 其 功率 通常 小 于 噪声 功率 但 大 于 射电 信和 号 功率 。 

通过 以 上 分 析 可 知 ， 如 何 抑制 射电 观测 中 弱 于 噪声 的 干扰 信号 是 本 文 讨 论 的 主要 问题 与 难点 。 对 
于 导航 信号 这 类 弱 于 噪声 的 干扰 信号 而 言 ， 完 全 获取 信和 号 的 先 验 信息 难度 较 大 ， 但 可 利用 信和 号 本 身 的 
村 点 来 抗 干扰 ， 如 空间 特性 和 相关 特性 。 同 时 ， 抗 干扰 方法 的 实施 难免 对 射电 信和 号 处 理 有 一 定 影响 ， 
因此 抗 干 扰 系统 的 建立 既 需 要 有 效 利用 干扰 信号 本 身 的 特点 ， 同 时 也 需 尽量 避免 对 射电 信和 号 处 理 的 影 
响 。 进 一 步 来 说 ， 有 效 地 利用 干扰 信号 的 空间 特性 需要 基于 多 通道 接收 的 方式 ， 文 [1-2] 利 用 干涉 
仪 和 相 控 阵 天 线 自 身 多 通道 的 特性 提取 干扰 信号 的 子 空间 ， 通 过 正 交 投影 的 方式 消除 干扰 信号 ， 然 而 
由 于 信号 空间 与 干扰 空间 不 垂直 ， 导 致 对 观测 信号 影响 较 大 。 为 了 利用 干扰 信号 的 相关 性 ， 借 助 辅 助 
天 线 是 常用 方式 。 文 [3-4] 通 过 辅助 天 线 估 计 更 为 精准 的 干扰 子 空 间 ， 并 运用 如 正 交 子 空间 、 互 子 空 
间 投 影 和 多 阵列 劳 澄 相 消 等 方式 抑制 干扰 信号 。 对 于 大 型 单 射电 天 线 而 言 ， 要 利用 干扰 信号 的 空间 特 
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性 和 相关 性 ， 则 需要 借助 多 个 辅助 天 线 。 文 [5] 从 雷达 信号 处 理 角度 出 发 ， 基 于 特征 子 空间 的 自 适应 
旁 狼 相 消 方法 ， 利 用 辅助 阵列 中 强 干扰 信和 号 的 特征 子 空 间 与 噪声 子 空 间 良 好 的 垂直 性 ， 消 除了 接收 信 
号 中 的 强 干扰 信号 。 对 于 抑制 射电 观测 中 弱 于 噪声 的 导航 信号 而 言 ， 需 根据 射电 信号 处 理 流程 重新 构 
建 抗 干扰 系统 框架 ,并 且 改 进 自 适应 旁 办 相 消 方法 有 效 抑制 导航 信和 号。 同时 考虑 到 射电 观测 数据 应 
用 ， 应 给 出 新 的 评价 指标 来 评估 干扰 抑制 程度 和 信号 损失 程度 。 此 外 ,实际 情况 中 射电 天 线 ( 主 天 
线 ) 口 径 通 常 较 大 ， 使 架设 辅助 阵列 离 主 天 线 较 远 ， 因 此 由 于 空间 位 置 产 生 的 信号 时 延 导 致 导航 信号 
相关 性 减少 ， 需 要 考虑 时 延 补偿 的 问题 。 

综 上 所 述 ， 本 文 首先 依据 射电 信号 处 理 流程 构建 了 基于 辅助 天 线 的 射电 抗 干扰 系 统 框架 ， 并 从 射 
电 观 测 数据 应 用 出 发 ， 给 出 了 干扰 抑制 度 和 信号 损失 度 两 个 指标 。 同 时 本 文 在 精确 补偿 主 天 线 接收 信 
号 时 延 后 ， 利 用 抑制 干扰 子 空间 蝶 声 分 量 的 自 适应 旁 办 相 消 方法 来 抑制 导航 信号 的 影响 。 最 后 利用 蒙 
特 卡 洛 仿 真 方式 ， 讨 论 了 不 同 主 天 线 和 辅助 阵列 间距 ( 即 主 、 辅 天 线 距 离 )、 干 扰 波 达 角 度 、 观 测 时 
间 、 辅 助阵 元 口径 和 阵 元 数量 对 抗 干扰 性 能 的 影响 。 


1 基于 辅助 阵列 的 射电 抗 干 扰 系 统 框 架构 建 


在 大 型 单 射电 天 线 中 ， 利 用 导航 信号 的 空间 特性 和 相关 性 抑制 干扰 ， 需 要 借助 多 通道 的 辅助 天 线 
实现 ， 即 辅助 阵列 。 自 适应 劳 辩 相 消 作为 常用 的 抗 干扰 方法 之 一 ， 本 质 是 利用 辅助 阵列 与 主 天 线 中 干 
扰 信 号 的 相关 性 和 干扰 的 空间 特性 获得 辅助 阵列 的 最 优 权 矢 量 ， 使 合成 的 天 线 方向 图 在 干扰 方向 上 产 
生 零 点 ， 达 到 干扰 抑制 的 目的 。 此 外 ， 为 了 降低 抗 干扰 系统 对 射电 信和 号 处 理 的 影响 ， 抗 干扰 系统 应 独 
立 于 射电 信号 处 理 流程 之 外 。 根 据 以 上 分 析 ， 建 立 基 于 辅助 阵列 的 射电 抗 干 扰 系 统 框架 ， 如 图 1。 
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数据 相关 处 理 软件 应 用 


图 1 基于 辅助 天 线 的 射电 抗 和 干扰 系统 框架 


Fig. 1 Framework of the radio anti-interference based on auxiliary antenna 
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通常 射电 接收 机 带宽 比 干 扰 信 号 带宽 宽 得 多 ， 例 如 上 海天 文 台 65 m 射电 天 线 442~2 400 MHzl9 ， 云 
南天 文 台 40 m 射电 天 线 300~800 MHz ， 因 此 需要 选取 包含 导航 信号 的 中 频 信 号 y(t+7) (如 图 1 中 
“1? 线 ) 为 抗 干 扰 系统 输入 的 原始 射电 接收 信号 。 值 得 注意 的 是 ， 混 频 器 1 是 补偿 主 天 线 接收 机 中 混 
频 器 本 振 频 率 与 辅助 天 线 接收 机 中 心 频率 的 差异 。 为 了 降低 对 采集 设备 的 性 能 要 求 ， 中 频 信 号 yarr) 
辅助 阵列 (x1(t) ,…, xi(t) ) 的 信号 通过 带 通 滤波 器 选取 存在 导航 信号 的 频带 信号 ， 在 经 过 混 频 器 2 
和 低 通 滤波 器 处 理 后 进入 多 路 数据 采集 系统 对 信号 进行 采集 。 而 由 于 空间 位 置 造 成 的 时 延 需 要 通过 
时 延 滤 波 器 实现 精确 补偿 ， 即 得 到 时 延 补偿 后 的 主 天 线 信 号 y(n)。 补 偿 后 的 y(n) 与 辅助 阵列 信号 
(x(n), +, xi(n) ) 通 过 自 适应 旁 闪 相 消 方法 ( 抗 干扰 算法 处 理 ) 得 到 y(n)。 将 Y(n) 进 行 数 模 转换 、 
上 变频 等 一 系列 逆向 处 理 得 到 y(1) ， 之 后 与 通过 带 阻 滤波 器 的 y(t+7) 信 号 相 加 得 到 抗 干扰 后 的 信号 
Y(t+7)。 以 上 抗 干 扰 过 程 中 需要 解决 几 个 关键 问题 : 

(1) 指标 构建 : 与 雷达 、 通 信和 领域 不 同 ， 射 电 天 文 抗 干 扰 中 干扰 信号 通常 小 于 噪声 信号 。 因 此 有 
必要 着 眼 于 射电 天 文 数据 应 用 提出 评价 指标 来 评估 干扰 消除 程度 和 信号 损失 程度 。 

(2) 时 延 补偿 : 由 于 射电 天 文 望远镜 口径 通常 较 大 ， 利 用 辅助 天 线 抗 干 扰 时 ， 其 架设 位 置 离 射 电 
天 线 ( 主 天 线 ) 位 置 较 远 ， 因 此 需要 精确 补偿 信和 号 时 延 ， 提 高 导航 信号 的 相关 性 。 

(3) Ai RA. 完全 已 知 所 有 导航 卫星 的 先 验 信息 难度 较 大 ， 但 其 运行 轨道 相对 固定 易于 
预测 ， 因 此 可 利用 导航 信号 良好 的 空间 特性 消除 影响 。 自 适应 旁 六 相 消 方法 中 基于 特征 子 空间 类 的 方 
法 可 以 在 干扰 较 强 时 ， 有 效 避 免 估 计 误 差 的 影响 ， 获 得 辅助 阵列 的 最 优 权 矢 量 。 然 而 导航 信号 功率 通 
常 小 于 噪声 功率 ， 因 此 需要 改进 此 类 方法 。 


oo 2 eT ERT SE AMES HY A E ON LT ST Hot PCT A 


his 2.1 评价 指标 

在 射电 抗 干扰 中 ， 干 扰 信 号 功率 通常 小 于 噪声 功率 ,传统 的 性 能 指标 如 相 消 比 等 ， 无 法 有 效 地 刻 
画 干扰 和 信号 分 量 的 变化 程度 。 为 解决 此 问题 ， 本 文 从 射电 观测 数据 应 用 出 发 ， 以 辅助 天 线 接收 信号 
以 及 先 验 的 射电 信号 为 标尺 ， 分 别 衡量 了 抗 干 扰 前 后 主 天 线 接收 信号 中 干扰 信号 与 射电 信号 的 相关 程 
度 ， 构 建 了 干扰 抑制 度 G 和 信号 损失 度 Gs 表征 干扰 抑制 程度 与 信号 损失 程度 ， 其 定义 分 别 为 : 


© FAMAE C: 
= _ max[R,,(7,) ] _max|E[y(t)«"(t-7,) |} 
© «max[Rs(7)] max} BLY(t)e* (+ 7) 1 
7 7 (1) 
_F(0) max|E[I(t)I*(t + ty -7,)]} _ FCO) 
F(0,) max E[IG)I (t +t) —7,)]} F)? 
同 理 ， 信 和 号 损失 度 G: 
max| Rs (T; F( 6. 
max[ Ra(71)] PCs) 车 


max[ Rs (7,) | E F(0,) 


HP, R(T), RoC) 表示 抗 干扰 前 、 后 辅助 天 线 接收 信号 x 与 主 天 线 接收 信号 yO) 的 互相 关 函 数 ，; 
Rs,(71) 、Rs;(7,) 表示 抗 干扰 前 、 后 参考 射电 信号 S 与 主 天 线 接收 信号 y(y) 的 互相 关 函 数 ;“ ~” 表 示 
抗 干扰 后 变化 的 量 ; max 表示 取 最 大 值 。 通 过 (1) 式 和 (2) 式 可 知 ， 干 扰 信 号 与 射电 信号 相关 性 的 变 
化 可 直接 反映 方向 图 在 干扰 和 信号 方向 上 的 变化 。 
2.2 ”时 延 补偿 

为 了 获取 良好 的 滤波 器 通 带 特性 ， 本 文采 用 基于 非 迭 代 的 加 权 最 小 二 乘 的 可 变 分 数 时 延 全 通 滤 波 
AES) BU GERMS)  ， 滤 波 器 结构 基于 farrow 结构 !"1 ， 该 结构 基于 不 同 准则 通 近 理想 的 滤波 器 获得 
滤波 右 系 数 后 ， 利 用 高 阶 分 数 时 延 的 多 项 式 近 似 得 到 滤波 器 系数 。 此 方式 即 可 将 分 数 时 延 与 滤波 器 系 
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数 分 离 ， 使 整个 滤波 右 的 传递 函数 简化 为 分 数 时 延 对 多 个 子 滤 波 絮 的 加 权 和 ， 因 此 该 结构 大 大 减少 了 
求解 滤波 器 系数 的 复杂 度 。 加 权 最 小 二 乘 滤波 方法 求解 滤波 器 系数 的 主要 思想 是 在 一 定 通 带 范围 和 时 
延 因 子 (分 数 时 延 与 采样 频率 的 乘积 ) 范围 内 ， 使 可 变 分 数 时 延 滤波 器 的 频率 响应 与 理想 的 分 数 时 延 
频率 响应 函数 误差 最 小 ， 即 


am 0.5 


minJ„(B)=min] f Wew, p)le,l@, p)|’dodp, (3) 


其 中 ,如 是 由 滤波 器 系数 组 成 的 矩阵 ; Woo, p) AERIAL; ey VAY ARN HE UE I tt EE 
分 数 时 延 频率 响应 的 误差 函数 ; p 为 除去 整数 部 分 的 分 数 时 延 因 子 。 不 难 发 现 (3) 式 中 存在 积分 过 程 ， 
利用 泰勒 级 数 展开 可 简化 积分 过 程 ， 即 通过 非 迭 代 的 方式 ， 快 速 获 得 滤波 器 系数 。 图 2(a) 和 图 2(b) 
分 别 是 分 数 时 延 因 子 p 变化 范围 为 -0. 5~0.5， 归 一 化 频率 o 变化 范围 为 0~0.97 时 ， 基 于 加 权 最 小 
二 乘 滤波 的 可 变 分 数 时 延 全 通 滤波 器 输出 的 相位 时 延 和 与 理想 的 分 数 时 延 频 率 响 应 函数 误差 结 


可 变 分 数 时 延 可 变 分 数 时 延 滤波 器 -频率 响应 误差 
x107 

tel 

EO z 

J = 
© 
ies 
| ee 
b 图 2 WLS 可 变 分 数 时 延 滤波 器 结果 
CO Fig.2 Results of the WLS-VFD filter 
N 
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确 的 分 数 时 延 。 


2.3 Ai DIA TA 

AERAN A A oC UE HR ERE, TRET E TA) AY FE MAT A E 
ræ 利用 干扰 子 空 间 与 噪声 子 空间 的 正 交 特性 ， 将 权 值 天 量 向 干扰 子 空间 投影 得 到 最 优 权 值 矢量 ， 可 见 准 
Q 确 佑 计 干 扰 子 空间 十 分 重要 。 辅 助阵 列 的 接收 信号 主要 是 干扰 信号 和 噪声 信号 ， 可 通过 佑 计 辅 助阵 列 
的 协 方差 矩阵 进而 获得 干扰 子 空 间 和 噪声 子 室 间 。 假 设 存在 p 个 互 不 相关 的 窄带 干扰 ， 协 方差 矩阵 R, 
特征 分 解 如 下 式 : 


R, =U,0,0' + UDD!, (4) 
Hop, R, Seth k(khEp) SABRE TORK EW N 的 接收 信号 估计 的 协 方差 矩阵 ;六 和 Dy 是 分 别 由 
前 bp 个 大 特征 值 和 后 k-p 个 小 特征 值 构成 的 对 角 阵 ;上 和 Uy 分 别 是 D, 和 Dy 对 应 特征 向 量 构成 的 矩 
阵 ， 其 列 矢量 分 别 张 成 的 空间 为 干扰 子 空 间 和 噪声 子 空间 。 

在 射电 天 文 观 测 中 消除 导航 信号 时 ， 导 航 信和 号 功率 通常 小 于 噪声 功率 ， 使 得 干扰 子 空间 受 噪 声 的 
扰动 较 大 ， 即 前 p 个 大 特征 值 和 后 -p 个 小 特征 值 差 异 较 小 ， 因 此 无 法 准确 地 获取 干扰 子 空间 。 为 此 
本 文 提出 改进 传统 的 基于 特征 子 空间 自 适 应 旁 闪 相 消 方法 ， 其 核心 思想 是 利用 噪声 的 平稳 特性 ， 通 过 
后 几 个 特征 值 估计 干扰 子 空间 中 的 噪声 分 量 并 剔除 ， 增 强 干 扰 子 空间 与 噪声 子 空 间 的 正 交 性 。 本 文 从 减 
小 噪声 对 干扰 子 空间 的 影响 和 消除 时 延 造成 导航 信号 相关 性 下 降 这 两 点 出 发 ， 估 计 最 优 权 值 矢 量 w,， 
的 步骤 为 : 

(1) 首 先 根据 导航 信号 入 射 方向 ( 即 波 达 角 ) 和 主 天 线 与 辅助 阵列 参考 阵 元 间 的 空间 位 置 估计 信号 
时 延 r; 
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(2) 根 据 信 号 时 延 7r 和 采样 频率 计算 时 延 因 子 1+p， 通 过 补 零 对 整数 时 延 了 补偿， 通过 2. 2 中 的 可 
变 分 数 时 延 滤 波 需 对 分 数 时 延 p 补偿 ; 
(3) 通 过 有 限 长 度 的 快 拍 数 估计 辅 助阵 列 的 协 方 差 矩 阵 和 主 天 线 与 辅助 阵列 的 相关 矢量 ， 即 
Ry 、R 
(AAT, 特征 分 解 如 (4) 式 ， 取 后 bcp 个 小 特征 的 平均 值 X=i， GRE TF 
T, R, 减 去 该 平均 值 构成 的 对 角 阵 ， 得 到 去 除 噪声 的 玉 : 
入 -| Yale (5) 


k TP i=p+1 


(5) 对 六 特征 分 解 ， 得 到 Ù, D, (6) SEE, B 
Won = UD U Ry è (6) 


通过 以 上 方式 可 将 干扰 子 空间 的 噪声 分 量 消除 一 部 分 ， 使 干扰 子 空间 与 噪声 子 空间 有 更 好 的 正 交 
性 ， 进 而 获得 更 优 的 抗 干扰 性 能 。 


3 仿真 结果 分 析 


3.1 参数 设置 

本 文 以 FAST 30 m 模型 脉冲 星 天 线 为 主 天 线 模型 "" 。 观 测 不 同 的 射电 信号 ， 其 信号 形式 有 较 大 
差异 ， 无 法 用 统一 的 形式 表示 射电 信号 ， 因 此 本 文 用 功率 谱 流 量 密度 为 1 mJy 的 随机 信号 代替 射电 信 
号 ; 辅助 天 线 是 由 无 方向 性 阵 元 构成 的 均匀 线 阵 ， 且 阵 元 间距 设置 为 0.09m。 以 调制 C/A 码 B3 频 点 
(1268. 52 MHz) 的 北斗 导航 信号 为 干扰 信号 。 

为 了 定量 地 对 比 时 延 补偿 对 本 文 方法 的 影响 和 本 方法 较 传统 的 基于 特征 子 空间 自 适 应 旁 办 相 消 方 
法 性 能 的 影响 ,分 别 对 比 了 在 时 延 补偿 和 时 延 未 补偿 情况 下 的 本 文 方法 与 基于 特征 子 空间 的 自 适 应 旁 
办 相 消 方法 ， 并 讨论 了 不 同 主 天 线 和 辅助 阵列 间距 ( 即 主 、 辅 天 线 距离 )、 干 扰 波 达 角 度 、 观 测 时 间 、 
辅助 阵 元 口径 和 阵 元 数量 对 抗 干 扰 性 能 的 影响 ， 参 数 设置 如 表 1。 此 外 ， 本 文 利用 蒙特 卡 洛 模 拟 方式 
消除 随机 误差 对 算法 评价 可 靠 性 的 影响 。 根 据 大 数 定理 可 知 ， 在 置信 和 度 w= 0.03， 方 差 为 0.074 3( 辅 
助 天 线 噪声 标准 差 )， 蒙 特 卡 洛 次 数 为 200 KIN, 蒙特 卡 洛 误差 为 a=0.015 8， 因 此 蒙特 卡 洛 次 数 为 
200 次 时 可 以 满足 比较 算法 性 能 的 要 求 。 

表 1 仿真 参数 对 比 
Table 1 Comparison of simulation parameters 


主 、 辅 天 线 距 离 /m 阵 元 直径 /m ”干扰 波 达 角度 (9)/。 ”观测 时 间 /s 辅助 阵 元 数量 


图 3(a)(b) 30 0.5 0~5 0. 01 10 
图 4(a)(b) 30 0.5 1.5 le-5~ 5e-3 10 
图 4(c)(d) 30 0.1~5 1.5 0. 01 10 
图 $(a)(b) 30 0.5 1.5 0.01 3~15 
HI 5(c) (d) 25~45 0.5 1.5 0. 01 10 


3.2 结果 分 析 
3.2.1 主 天 线 与 辅助 阵列 距离 的 影响 

主 天 线 与 辅助 阵列 中 参考 阵 元 的 距离 远近 导致 接收 信号 的 时 延 不 同 进而 影响 导航 信号 的 相关 性 。 
因此 本 文 讨论 了 主 、 辅 天 线 距离 对 方法 性 能 的 影响 ， 如 图 3。 

结合 表 1 和 图 3(a) 可 知 ， 随 着 主 天 线 与 辅助 阵列 距离 的 增加 ， 时 延 补偿 的 本 文 方法 干扰 抑制 度 比 
时 延 补偿 和 未 补偿 的 自 适 应 旁 澄 相 消 方法 高 15 dB 左右 ， 同 时 ， 时 延 未 补偿 的 本 文 方法 干扰 抑制 度 ， 由 
于 时 延 的 影响 随 主 、 辅 天 线 距离 增加 而 下 降 ; 另外 ， 无 论 时 延 是 否 补偿 ， 本 文 方法 和 基于 特征 子 空间 的 
自 适 应 旁 注 相 消 方法 ， 随 着 主 天 线 与 辅助 阵列 距离 增加 ， 信 号 损失 度 相差 不 大 , 约 1dB， 如 图 3(b)。 
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图 3 不 同 主 天 线 与 辅助 阵列 距离 的 方法 性 能 对 比 


Fig. 3 Performance comparison with different distance between the main-antenna and the auxiliary antennas as the variable 


3.2.2 波 达 方 向 和 阵 元 口径 的 影响 
不 同 来 向 的 导航 信号 由 于 方向 图 的 影响 ， 功 率 差 异 很 大 ， 导 致 干扰 信号 (导航 信和 号) 在 主 天 线 中 
干 品 比 差异 也 很 大 ; 另 一 方面 ， 阵 元 口径 影响 干扰 信号 的 接收 功率 。 因 此 本 文 讨论 了 不 同 信号 波 达 角 
0( 与 主办 轴 的 夹 角 ) 和 阵 元 口径 对 方法 性 能 的 影响 ， 如 图 4。 
干扰 抑制 度 /Gi 号 损失 度 (Gs 
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图 4 不 同 波 达 角 和 辅助 阵 元 口径 的 方法 性 能 对 比 


Fig.4 Performance comparison with direction of arrival as the variable and aperture of the auxiliary antennas as the variable 
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结合 表 1 和 图 4(a) 可 得 ， 随 着 导航 信号 波 达 角 远 离 主 欠 ， 其 信号 功率 也 随 方向 图 变化 而 减 小 。 
时 延 补偿 和 时 延 未 补偿 的 本 文 方法 干扰 抑制 度 也 随 之 减 小 ， 但 是 时 延 补偿 的 本 文 方法 比 自 适应 旁 泊 相 
消 方法 干扰 抑制 度 高 7~20 dB， 而 时 延 未 补偿 的 本 文 方法 仅 在 信号 波 达 角 小 于 2.5° 时 比 自 适应 劳 准 相 
消 方法 干扰 抑制 度 高 ;图 4(b) 中 导航 信号 波 达 角 在 1. 5 范围 内 信号 损失 度 比 劳 注 方向 上 更 大 ， 是 因 
为 合成 的 方向 图 使 主流 零 陷 ， 零 陷 程 度 越 深 ,信号 的 损失 越 大 ， 因 此 无 论 时 延 是 否 补偿 ， 本 文 方法 信 
号 损失 度 比 自 适应 旁 准 相 消 方法 高 约 3dB。 当 导航 信号 从 旁 淮 甚至 远 旁 办 进入 接收 机 时 ， 旁 办 出 现 零 
陷 ， 因 此 信和 号 损失 度 Cs 逐渐 趋 于 0。 

结合 表 1 和 图 4(c) 可 得 ， 随 着 辅助 阵 元 口径 增加 (等 价 于 干 噪 比 的 增加 ) ， 时 延 补偿 的 本 文 方法 
和 自 适应 旁 辩 相 消 方法 干扰 抑制 度 G 均 不 断 提高 ， 由 于 口径 的 增 大 使 得 干扰 信号 接收 功率 增 大 。 时 
延 补偿 的 自 适应 旁 兴 相 消 方法 在 辅助 阵 元 口径 大 于 1.4 m 时 , 干扰 抑制 度 超过 时 延 未 补偿 的 本 文 方 
法 ,但 其 干扰 抑制 度 比 时 延 补偿 的 本 文 方法 的 G1 低 7~10 dB。 随 着 阵 元 口径 的 增加 ， 本 文 方法 与 自 
适应 旁 兴 相 消 方法 的 信号 损失 度 相差 2 dB 内 ， 如 图 4(d) 。 
3.2.3 观测 时 间 和 阵 元 数量 的 影响 

观测 时 间 是 指 参与 估计 RR,、RR,, 运 算 的 快 拍 数量 ( 即 接收 信号 数据 长 度 )N 所 对 应 的 时 间 ， 不 同 的 
观测 时 间 影 响 相关 和 矩阵 的 佑 计 精 度 ， 而 阵 元 数量 影响 合成 方向 图 在 干扰 来 向 零点 的 零 陷 程度 和 可 消除 
干扰 个 数 ， 因 此 本 文 讨论 了 不 同 观测 时 间 和 阵 元 数量 对 方法 性 能 的 影响 ， 如 图 5。 

结合 表 1 和 图 5(a) 可 得 ， 随 着 观测 时 间 的 增加 ， 时 延 补偿 的 本 文 方法 干扰 抑制 性 能 不 断 提高 ， 
干扰 抑制 度 比 自 适应 旁 兴 相 消 方 法 高 10~18 dB， 且 比 时 延 未 补偿 的 本 文 方法 干扰 抑制 度 高 10~12 dB, 
相 比 自 适应 旁 泊 相 消 方法 ， 时 延 补偿 的 本 文 方法 的 信号 损失 度 Cs 高 约 1dB， 如 图 5(b)。 
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图 5 不 同 观测 时 间 和 阵 元 数量 的 方法 性 能 对 比 


Fig.5 Performance comparison between different observation time and between different numbers of auxiliary antenna 
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结合 表 1 和 图 5(c) 可 知 ， 随 着 阵 元 数 的 增加 ， 两 种 抗 干 扰 方 法 的 干扰 抑制 度 都 会 随 之 增加 ,但 
时 延 补偿 的 本 文 方法 比 自 适应 旁 办 相 消 方法 干扰 抑制 度 C 高 约 14 dB， 且 比 时 延 未 补偿 的 本 方法 干扰 
抑制 度 高 约 4 dB。 此 外 ， 无 论 时 延 是 否 补偿 ， 本 文 方法 的 信号 损失 度 C 整体 比 自 适 应 劳 办 相 消 方法 
高 约 1dB， 如 图 5(d)。 
4 结 论 

本 文 针 对 导航 信号 干扰 L 波段 射电 观测 这 一 实际 情况 ,构建 了 基于 辅助 阵列 的 抗 干扰 系统 框架 ， 
并 从 射电 数据 应 用 出 发 ， 提 出 了 干扰 抑制 度 和 信和 号 损失 度 两 个 指标 。 利 用 导航 信号 良好 的 空间 特性 和 
相关 特性 ,改进 了 传统 的 基于 特征 子 空间 的 自 适 应 旁 辩 相 消 方法 ， 使 其 在 干扰 功率 小 于 噪声 功率 时 
依然 可 有 效 地 消除 干扰 。 此 外 ， 本 文 考虑 到 由 于 射电 天 线 与 辅助 阵列 空间 相对 位 置 造成 的 时 延 导致 
导航 信号 相关 性 降低 问题 。 利 用 基于 非 迭 代 的 加 权 最 小 二 乘 的 可 变 分 数 时 延 全 通 滤 波 器 实现 了 时 
延 的 精确 补偿 ， 提 高 了 导航 信和 号 的 相关 性 。 最 后 基于 所 提 指 标 ， 通 过 蒙特 卡 洛 仿真 方式 ， 对 比 了 本 文 
方法 与 自 适应 劳 闪 相 消 方法 ， 结 果 表 明 本 文 方法 的 干扰 抑制 度 提 高 了 10~20 dB， 但 信号 损失 度 增加 
约 1 dB, 
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Abstract: When navigation signals, as radio interferences, affect observation in L-band, they are 
different from other types of interferences in communication field or radar. Using the correlation interference , 
traditional methods are ineffective and they lack effective evaluation indexes, too, because the power of 
navigation signal is normally weaker than noise power. To solve the above problems, this paper constructs a 
system framework of radio astronomy interference based on radio signal processing, and raises two indexes; 
interference suppression and signal loss based on the application of the observation data. After accurate delay 
compensation, this paper proposes a method to eliminate the effects of the navigation signals, using the 
adaptive side-lobe cancellation that suppresses the noise component of interference subspace. From the results 
of the simulation, conclusion can be drawn that the proposed method can effectively suppress the navigation 
signal in the observations, and the performance is better than that of traditional adaptive side-lobe cancellation 
based in feature subspace. 

Key words; Radio astronomy interference; Adaptive side-lobe cancellation; Navigation signal; Variable 


fractional delay filter 


